Conductance of a quantum point contact based on spin-density-functional theory 
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We present full quantum mechanical conductance calculations of a quantum point contact (QPC) 
performed in the framework of the density functional theory (DFT) in the local spin-density ap- 
proximation (LDA). We start from a lithographical layout of the device and the whole structure, 
including semi-infinitive leads, is treated on the same footing (i.e. the electron-electron interac- 
tion is accounted for in both the leads as well as in the central device region). We show that a 
spin-degeneracy of the conductance channels is lifted and the total conductance exhibits a broad 
plateau-like feature at ~ 0.5 x 2e 2 /h. The lifting of the spin-degeneracy is a generic feature of 
all studied QPC structures (both very short and very long ones; with the lengths in the range 
40 < I < 500 nm). The calculated conductance also shows a hysteresis for forward- and backward 
sweeps of the gate voltage. These features in the conductance can be traced to the formation of 
weakly coupled quasi-bound states (magnetic impurities) inside the QPC (also predicted in previous 
DFT-based studies). A comparison of obtained results with the experimental data shows however, 
that while the spin-DFT based "first-principle" calculations exhibits the spin polarization in the 
QPC, the calculated conductance clearly does not reproduce the 0.7 anomaly observed in almost all 
QPCs of various geometries. We critically examine major features of the standard DFT-based ap- 
proach to the conductance calculations and argue that its inability to reproduce the 0.7 anomaly 
might be related to the infamous derivative discontinuity problem of the DFT leading to spurious 
self-interaction errors not corrected in the standard LDA. Our results indicate that the formation of 
the magnetic impurities in the QPC might be an artefact of the LDA when localization of charge is 
expected to occur. We thus argue that an accurate description of the QPC structure would require 
approaches that go beyond the standard DFT+LDA schemes. 



I. INTRODUCTION 

Experimental evidence of an additional 0.7 x 
2e 2 /h conductance feature in quantum point contacts 
(QPCs)L2 iMi i L £ iL £i2il£ has generated enormous theoret- 
ical activity during recent decade&i&ii^^^^^^^. 
While a usual conductance quantization in terms of 
2e 2 /h can be successfully explained in a one-electron 
picture 1 ^, no consensus on the origin of the 0.7 anomaly 
has been reached so far. Experimental data, based 
on magnetic field dependence of the plateau position 1,2 
and observation of the zero bias anomaly^ clearly 
point out at spin origin of this effect. The spin 
origin was adopted in many theories which include, 
just to name a few of them, a spontaneous spin- 
polarization inside the constriction of the QPC origi- 
nated from the exchange-correlation interactio n 11 ' 12 ! 13 ' 14 , 
the ferromagnetic-antifcrromagnetic exchange interac- 
tion with a large localized spin™, a formation of polar- 
ized quasi-bound states and Kondo effec t 15 ' 16 , as well as 
several others. However, none of the existing theories re- 
produce quantitatively the 0.7 conductance anomaly ob- 
served in almost all QPC^SA^L^.. Note that some 
of approaches based on phenomenological models do re- 
cover the 0.7 feature in the conductance^^. Such the 
approaches, while providing an important insight for an 
interpretation of the experiment, are not, however, able 
to uncover the microscopic origin on the observed effect. 

A powerful technique capable of providing detailed in- 
formation about electronic and transport properties on 
the microscopic level is based on the " first principle" den- 



sity functional theory (DFT) approach. The spin DFT 
calculations addressing the 0.7 anomaly in the QPC have 
been presented by various group o 11 ' 12 ' 13 ' 14 ' 15 ' 17 . How- 
ever, many of these calculations show rather conflicting 
results. For example. Refs. 1 5 1 (i. \ l\ attribute the 0.7 
anomaly to formation of the localized magnetic moment 
in the QPC, whereas such quasi-bound states are not re- 
covered in Refs. mum!. In contrast, Ref. relates 

the 0.7 anomaly to multiple metastable spin-polarized 
solutions. The recent study of Jaksch et al~ suggests 
that the spin polarization is absent in short QPCs and 
is increased with the increase of the length of the QCP. 
This seems to contradict the results reported by Rejec 
and Meir— who find that a strong spin polarization due 
to the magnetic impurity formation is a generic feature 
seen in short as well as in long QPCs. 

Motivated by the above mentioned conflicting results 
and conclusions, in this paper we perform "first princi- 
ple" full quantum mechanical transport calculations of 
the conductance of the QPC based on the spin DFT in 
the local spin-density approximation (LDA). [We use a 
notation "first principle" in a sense that we start from a 
lithographical geometry of the gates and a heterostruc- 
ture layout]. The main features of our calculations and 
the relation of our model to the DFT-based calculations 
reported previously are as follow. In the present paper 
the QPC is viewed as an inherently open system where 
the solution of the Schrodinger equation consists of con- 
tinuous scattering states. This is in contrast to Refs. 
[ITIIT3I treating a QPC as a closed system where the so- 
lution of the Schrodinger equation is represented by a 
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discrete set of the eigenf unctions. In our calculations we 
start from a lithographical layout of the device, which 
is in contrast to Ref. using a simplified model of an 
external confinement potential. (Note that importance 
of an accurate treatment of the confinement potential is 
stressed in Ref. [H where the phenomenology and ex- 
perimental data indicate that the 0.7 feature depends 
strongly on the potential profile of the contact region). 
Our model (but not the computational technique) is con- 
ceptually similar to the one used by Rejec and Meir— (see 
also Ref. [TH reporting similar findings) where a QPC is 
treated as an open structure and where a realistic model 
for the external confinement potential due to metallic 
gates is utilized. In their study Rejec and Meir focus on 
the calculation of the electron density, predicting the for- 
mation of a localized spin-degenerate quasi-bound state 
(magnetic impurity). They did not, however, perform 
transport conductance calculations so the central ques- 
tion whether the 0.7-anomaly is indeed related to the for- 
mation of a magnetic impurity has remained unanswered. 
In the present paper we hence focus on the calculation 
of the spin-resolved conductance of the quantum point 
contact. 



The details of our model and the Hamiltonian are given 
in Sec. II. In Sec. Ill we present our method for the cal- 
culation of the electron density and conductance which 
is based on the self-consistent Greens function technique 
combined with the spin DFT in the local spin density 
approximation. This method is a generalization of a 
method of Ref. for the case of spin. The results are 
presented in Sec. IV. Our calculations reconfirm the find- 
ing reported by Rejec and Meir— concerning the forma- 
tion of the localized spin-degenerate quasi-bound states 
(magnetic impurities) within the DFT+LDA approach. 
However, the calculated conductance clearly does not re- 
produce the 0. 7 anomaly observed in almost all QPCs of 
various geometries. Instead, the total conductance shows 
a broad feature peaked at 0.5 x 2e 2 jh. A similar feature 
is also present in the range of the gate voltages where 
a second step in the conductance develops. The calcu- 
lated conductance also shows a hysteresis for forward- 
and backward sweeps of the magnetic fields. (We stress 
that all these results are generic; we studied QPS with 
lengths in the range 40 - 500 nm and electron densities 
in the leads in the range I0 15 m -2 — 4 • 10 15 m~ 2 , with 
very similar results). In Sec. V we discuss the obtained 
results and critically examine the DFT+LDA-based con- 
ductance calculations. We suggest that the failure of the 
DFT+LDA approach to reproduce quantitatively the 0.7 
anomaly may be due to the lack of the derivative dis- 
continuity in the standard LDA approximation (leading 
to uncorrected self-interaction errors) for the case when 
localization of charge is expected to occur, so that the 
magnetic impurity formation may be an artefact of the 
DFT+LDA approach due to the spurious self-interaction. 
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FIG. 1: (color online) (a) Structure of a quantum point con- 
tact typically fabricated by a split-gate technique. An inter- 
nal computational region is attached to two semi-infinite leads 
which serve as electron reservoirs. The gates of the quantum 
point contact form a constriction of the length I and the width 
w; Wi e ad is the distance between the gates defining the leads, 
(b) A representative calculated electron density. 



II. MODEL 

We consider a quantum point contact (QPC) placed 
between two semi-infinite electron reservoirs. A 
schematic layout of the device is illustrated in Fig. 
[lja). Charge carriers originating from a fully ion- 
ized donor layer form the two-dimensional electron gas 
(2DEG) which is buried inside a substrate at the 
GaAs/ALGai-^As heterointerface at the distance b from 
the surface. Metallic gates are situated on the top of the 
heterostructure and define the quantum point contact as 
well as electron reservoirs that are are represented by 
uniform quantum wires of the infinite length. 

The Hamiltonian of the whole system (QPC + leads) 
within the Kohn-Sham formalism can be written in the 
form H = £ CT i? CT , 



H" = 



2m 



-V 2 + V a (r), 



(1) 



r = (x, y), m* = 0.067to c is the GaAs effective mass and 
a stands for spin-up, t, and spin-down, j, electrons. The 
first term in ([T]) is the kinetic energy of an electron while 
V a {r) = V con f(r) + V H (r) + V° c (r) is the total confin- 
ing potential which is a sum of the electrostatic confine- 
ment potential, the Hartree potential, and the exchange- 
correlation potential, respectively. The electrostatic con- 
finement V conf (r) = Vgates(r) + V d + V S chottky in- 
cludes contributions from the top gates, the donor layer 
and the Schottky barrier. The explicit expressions for the 
potentials V ga tes (r) and Vdonors are given in Refs. [2l| and 
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l2a the Schottky barrier is chosen to be Vschottky = 0.8 
eV. The Hartree potential is written in a standard form 



Mr) 



4ir£ e r 



dr 'n(r') 



1 



1 



|r-r'| ^Ir-r'l 2 +46 2 

where n(r) = J2a n(T ( r ) is the total electron density and 
the second term describes the mirror charges placed at 
the distance b from the surface, e r = 12.9 is the dielectric 
constant of GaAs, and the integration is performed over 
the whole device area including the semi-infinite leads. 

The exchange-correlation potential V£ c (r) in the lo- 
cal spin density approximation is given by the functional 
derivative^ 



(3) 



For e xc we have employed two commonly used parame- 
terizations, by Tanatar and Ceperley2i and by Attacalite 
et al2^. These two parameterizations give very similar re- 
sults. All the results presented below correspond to the 
parameterization of Tanatar and Ceperlej«2£. 



III. METHOD 

The central quantity in transport calculations is the 
conductance. In the linear response regime, it is given 
by the Landauer formula G — J2 a G a r^- 



G a = -- 



dE T a (E) 



Of (E - E F ) 
dE 



(4) 



where T a (E) is the total transmission coefficient for the 
spin channel a, f(E — Ep) is the Fermi-Dirac distribution 
function and Ep is the Fermi energy. In order to calculate 
T a (E) we use the method developed in Ref. l27| for the 
spinless electrons and generalize it here for the presence 
of two spin channels. For the case of completeness, the 
main steps of our calculations are presented below. 

We discretize Eq. ([Tj and introduce the tight-binding 
Hamiltonian with the lattice constant of a — 4 nm. (Such 
a small a ensures that the tight-binding Hamiltonian is 
equivalent to the continuous Schrodinger equation.) The 
retarded Green's function is introduced in a standard 



way 26 , 



g a = (E - H a + irj)~ 



(5) 



The Green's function in the real space representation, 
Q" (r, r, E) , provides an information about the electron 
density at the site r,— 



r (r) = -^3 J dEg°(r,r,E)f(E-E F ). 



(6) 



Note that Q a (r, r, E) is a rapidly varying function of en- 
ergy. As a result, a direct integration along the real axis 



in Eq. ([6]) is rather ineffective as its numerical accu- 
racy is not sufficient to achieve a convergence of the self- 
consistent electron density. Because of this, we transform 
the integration contour into the complex plane £s[E] > 0, 
where the Green's function is much more smoother (see 
Ref. M for details). 

In order to calculate the Green's function of the whole 
system (QPC + leads) we divide it into three parts, the 
internal region and two semi-infinite leads, as shown in 
Fig. [JJa). The internal region consists of the QPC as well 
as a part of the leads. In order to link the internal region 
and the leads together the self-consistent charge density 
(and the potential) at both sides must be the same. To 
fulfill this requirement, we place the semi-infinite leads 
far away from the central region so that the gates defining 
the QPC do not affect the electron density distribution 
in the leads. (The distribution is determined solely by 
the lead gate voltage, see Fig. [TJ) As a result, it is 
fully justified to approximate the semi-infinitive leads by 
a uniform quantum wire. The self-consistent solution for 
the latter can be found by the technique developed in Ref. 
[28l . Eventually, the Green's function of the whole system 
is calculated by linking the surface Green's function for 
the semi-infinite leads (calculated using the technique of 
Ref. HH) and the Green's function of the internal region 
with the help of the the Dyson equation^. 

All the calculations described above are performed self- 
consistently in an iterative way until a converged solution 
for the electron density and potential (and hence for the 
total Green's function) is obtained. Having calculated 
the total self-consistent Green's function, the scattering 
problem is solved where the scattering states in the leads 
(both propagating and evanescent) are obtained using 
the Green's function technique of Ref. [28|. 

Having calculated the Green's function we can find the 
local density of states (LDOS) in the QP(M 



LINOS' 7 {x,E) = -ig? J dy g a (x,y,E). 



(7) 



The self-consistent solution in quantum transport or 
electronic structure calculations is often found using a 
"simple mixing" method. It is the most robust and reli- 
able algorithm with only one disadvantage, namely a low 
convergence rate. The charge density on the ( i + l)-th 
iteration loop is updated through the input n\ n (r) and 
output n° ut (r) densities on the previous i-th iteration 



^ 1 (r) = (l-e)nf(r) + enr t (r) ! 



(8) 



with e being a small constant ~ 0.1 — 0.01. It is typically 
needed ~ 200 — 2000 iteration steps to achieve our con- 
vergence criterium for the relative density update on the 
z-th iteration step, 



< 10" 



(9) 



where n — J n(r) dr is the total electron number in the 
internal region during the i-th iteration step. 
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In order to trigger a spin-polarized solution we apply at 
early iterations a small parallel magnetic field, B = 0.05 
T. The final solution is driven mostly by the exchange 
interaction which is about four orders of magnitude larger 
that the Zeeman energy at B = 0.05 T. 

The self-consistent calculation are performed for dif- 
ferent gate voltages. To facilitate the calculations we use 
a solution from the previous value of the gate voltage as 
an initial guess for the subsequent one. It is worth not- 
ing that the modified Broyden's second method^, which 
can greatly reduce the number of iterations for the case of 
spinless electrons 2 ^, does not lead to reliable convergent 
results for the case of the spin degree of freedom. 



IV. RESULTS 

We calculate the conductance of a split-gate QPC with 
following parameters representative for a typical experi- 
mental structure. The 2DEG is buried at b = 70 nm be- 
low the surface (the widths of the cap, donor and spacer 
layers are 24 nm, 36 nm and 10 nm respectively), the 
donor concentration is 0.6 • 10 24 m~ 3 . The width of the 
semi-infinitive leads is wi eac i = 400 nm, and the width of 
the QPC constriction is w = 100 nm, see Fig. [T](a). The 
width w is kept constant throughout the paper, while 
the length is varied in the range I = 40 — 500 nm. The 
gate voltage applied to the gates defining the leads is 
Viead = —0.73 V. With these parameters of the device 
there are 18 spin-up and 18 spin-down channels available 
for propagation in the leads and the electron density in 
the center of the leads is n] ead = n\ ead = 2.2 • 10 15 m -2 . 
Such a large number of the channels in the leads makes 
its DOS to mimic a smooth DOS of the two-dimensional 
electron gas. The temperature is fixed at T = 0.2 K for 
all results presented below. 

Figure [2] shows the conductance of QPCs of differ- 
ent constriction lengths I ranging from I = 50 nm (a 
very short QPC) to / = 500 nm (a long quantum wire- 
type QPC). The conductance of all QPCs shows a broad 
plateau-like feature at 0.5 x 2e 2 /h. As the length of the 
constriction increases, a dip following 0.5-plateau starts 
to develop in the QPC conductance. An inspection of 
the spin-resolved conductance demonstrates that 0.5 fea- 
ture corresponds to the transmission of only one spin 
channel (say, spin- up), whereas the second (spin-down) 
conductance channel is totally suppressed. For long con- 
strictions (I > 300 nm) the 0.5 plateau starts to "wear 
down" transforming into a broad feature whose maximal 
amplitude is less than 0.5 x 2e 2 /h. If the constriction 
is sufficiently long (I > 100 nm), a conductance plateau 
at ~ 1.5 x 2e 2 /h starts to develop and a conductance 
dip following the 1.5- plateau also starts to emerge as the 
length of the constriction increases. 

To shed light on a microscopic origin of the 0.5 x 2e 2 /h 
conductance feature and the suppression of the spin- 
down channel let us inspect the charge density and energy 
structure of the QPC. Figure [3] shows the charge densities 
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FIG. 2: (color online) Conductance of the quantum point 
contact of different lengths 50 < I < 500 nm as a function of 
the gate voltage V g . The geometrical width of the constriction 
is i« = 100 nm; the geometrical length I is indicated in the 
figures. The dashed line corresponds to the spin-unpolarized 
solution. 



and the LDOS in the QPC of the length of I = 200 nm 
for the case of one transmitted spin-up, G^ = e 2 /h, and 
totally blocked spin-down, G- 1 = 0, channel. Near the 
constriction entrance, the spin-polarized charge droplets 
(marked by dotted circles in Fig. 13(a)) are clearly visible. 
The spin polarization is caused by the exchange interac- 
tion which dominates the kinetic energy for low densi- 
ties. The integration of the electron density gives about 
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FIG. 3: (color online) Formation of quasi-bound states in 
the quantum point contact, (a) The charge density and (b) 
the local density of states are shown for the regime of one 
transmitted spin-up and totally blocked spin-down channel. 
The left and right columns correspond to spin-up and spin 
down electrons. The white dashed lines in (b) indicate the 
self-consistent Kohn-Sham potential along the cross section 
y — 0, V' J (x,y = 0). The charge droplets with approximately 
one electron trapped are denoted by dotted lines in (a). The 
geometrical length and width of the QPC are I — 200 nm and 
w = 100 nm, respectively; V g = —0.98 V. 



one particle in each droplet. The quasi-bound states for 
corresponding droplets can also be traced in the LDOS 
shown in Fig. [3fb). An inspection of the correspond- 
ing potential profile for the spin-up electrons reveals that 
the low- and high energy droplets correspond to respec- 
tively first and second quasibound states trapped in the 
double-barrier potential. At the same time the spin-down 
droplets are spatially separated by a distance ~ 100 nm. 
This spatial separation is also reflected in the shape of the 
total confining potential for the spin-down electrons that 
forms a wide tunnelling barrier, see Fig. [3] (b). Because 
of this barrier, the transmission probability for the spin- 
down conductance is negligibly small (G^ ~ 10~ 5 xe 2 /h). 

The quasi-bound states are not spatially fixed but 
gradually move during the sweep of the gate voltage, V g . 
Already in the pinch-off regime, two quasi-bound states 
are developed at both sides of the constriction. When 
V g becomes less negative they move towards each other 
and eventually merge. Because of the exchange inter- 
action, this occurs first for the spin-up state and then 
for the spin-down. As a result, the total conductance G 
shows quantization in units of e 2 //i, see Fig. [2] We stress 
that all the results reported above are generic; we studied 
QPCs with lengths in the range 40 - 500 nm and electron 
densities in the leads in the range 10 15 m~ 2 — 4 • 10 15 m~ 2 , 
with very similar results. We also stress that our calcu- 
lations with two parameterizations of Refs. [24ll25l for the 
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FIG. 4: (color online) Convergence of self-consistent calcu- 
lations for (a) a short QPC (I — 100 nm) and (b) a long 
QPC (I = 300 nm) (y-axis corresponds to the relative density 
update defined by Eq. (0). Solid and dashed lines denote 
spin-up and spin-down solutions, respectively. Insets show 
the charge densities at the constriction at different stages of 
iteration procedure. Note, that we used the convergence crite- 
ria for the relative density update < 10 -7 for all calculations 
presented in this paper. The geometrical width of the con- 
striction w — 100 nm; V g = —1.4 V. 



correlation potential give almost the same results. This 
is not surprising, because for the system at hand the cor- 
relation potential is an order of magnitude lower than the 
exchange potential. 

Our conclusions concerning formation of the localized 
quasi-bound states in the QPC agree well with earlier 
findings of Hirose et al^ and Rejec and Meir— , but do 
not support the conclusion of Jaksch et al^ that the spin 
polarization is absent in short QPCs and is increased with 
the increase of the length of the QCP. 

We now want to focus on an aspects of the calcula- 
tions concerning a choice of the criteria for termination 
of the iteration procedure, Eq. ([9]). Even though this 
aspect has a rather technical character, we feel that it 
is important to address it in detail, as the proper choice 
of the termination criteria is essential for reaching of the 
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FIG. 5: (color online) Conductance hysteresis in the quan- 
tum point contact for forward- and backward sweeps of the 
gate voltage (solid and dashed lines respectively). The dot- 
ted line corresponds to the unpolarized solution exhibiting no 
hysteresis behavior. Insets show the charge densities in the 
constriction at V a = —0.975 V. The parameters of QPC are 
I = 200 nm and w — 100 nm. 



spin polarized solution. Figure Q] shows a relative den- 
sity update as a function of the number of iterations for 
a relatively short (I — 100 nm) and a relatively long 
(I = 300 nm) constrictions. For the long constriction the 
spin polarized solution is obtained already during initial 
iteration steps and its character does not change with 
increase of the number of iterations, see Fig. [4] (b). On 
the contrary, for a short constriction, the spin polarized 
solution is obtained only during later stages of iterations, 
as no spin-polarized solution is present during the initial 
steps even when the relative density update decreases to 
~ 10~ 4 , see 0] (a). Note, that the spin-polarized solu- 
tion is energetically favorable and more stable than the 
spin-unpolarized one for all QPC studied. 

The spin-DFT calculations reported above show the 
presence of the spin-resolved quasi-bound states confined 
in an effective double-barrier potential. It is well known 
that the quasi-bound states often lead to the hysteretic 
behavior of the system. For example, the hysteresis is 
commonly observed in the I- V characteristics of the res- 
onant double-barrier tunnelling structures for forward- 
and backward voltage sweep o 31 i 32 . Such hysteretic be- 
havior is related to the presence of the quasi-bound 
state between the barriers and was reproduced in nu- 
merous calculations 32 , including those based on the DFT 
approach 33 . Our calculations demonstrate that similar 
hysteretic behavior is present in the conductance of the 
QPC. When the gate voltage is swept forward and back- 
wards, the conductance shows the prominent hysteresis 
in the transition regions between the plateaus, see Fig. 
[5j In the region of hysteresis, the system, depending on 
the history, can be in one of two different ground states 



as illustrated in the inset to Fig. [5l Because the hystere- 
sis is present only in the transition regions and is com- 
pletely absent in the spin-unpolarized solution, its origin 
is due to the exchange interaction. We would like to 
stress that existence of two distinct ground states as well 
as the hysteresis behavior of the system at hand does not 
contradict to the DFT approach as such. Owing to the 
nonlinearity of the problem, there might be more than 
one solution within the DFT framework 34 . Besides the 
double-barrier resonant tunnelling structures mentioned 
above 33 , the DFT approach was used for the description 
of the hysteretic behavior of quantum Hall states in the 
double-layered quantum well structures^., as well as hys- 
teresis and spin phase transitions in quantum wires in 
the integer quantum Hall regime^. We however, con- 
clude this section by noticing that we are not aware of 
any experimental reports on the observation of the hys- 
teresis phenomena in the QPC structures. 



V. DISCUSSION 

We start our discussion from a comparison of the spin 
DFT-based conductance calculations reported in the pre- 
ceding section with the experimental conductance of the 
QPC structures. The calculations show a pronounced 
plateau-like feature at ~ 0.5 x 2e 2 /h for all QPC struc- 
tures studied (both short and long ones). For longer 
QPCs (I > 150 nm), a dip following 0.5-plateau starts 
to develop in the conductance. At the same time, a 
new conductance plateau at ~ 1.5 x 2e 2 /h starts to de- 
velop and a conductance dip following 1.5-plateau also 
starts to emerge as the length of the constriction in- 
creases (see Fig. [5]). Our calculations also show a pro- 
nounced hysteretic behavior in the transition regions be- 
tween the plateaus for forward and backward sweeps of 
the gate voltage, Fig. \5\ On the contrary, the exper- 
imental data clearly show an anomaly in the conduc- 
tance around 0.7 x 2e 2 /h. [It should be stressed that the 
calculated 0.5 plateau corresponds to the complete sup- 
pression of one spin channel, whereas experimentally ob- 
served 0.7 anomaly means that both spin channels have 
to be present in the conductance]. Some experiments 8 
also indicate a feature around 1.7 x 2e 2 /h. As far as 
the hysteresis behavior is concerned, we are not aware 
of any experimental reports of this effect. The above 
comparison demonstrates that while the spin-DFT based 
"first-principle " calculations predict the spin polarization 
in the QPC structure, the calculated conductance clearly 
does not reproduce the 0. 7 anomaly observed in almost 
all QPCs of various geometries. 

In order to understand why the calculated conduc- 
tance fails to reproduce the 0.7 anomaly, let us critically 
analyze the major features of the DFT-based conduc- 
tance calculations. Our computation scheme relies on 
an approach that during recent years became a stan- 
dard tool for transport calculations in open electronic 
systems such as molecules, metallic wires and meso- 
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scopic conductor s 37 ' 38 ' 39 ' 40 ' 41 ' 42 . Its starting point is the 
Landauer-type formula where the conductance is calcu- 
lated using the Green's function or similar scattering 
techniques combined with the mean-field description of 
the electronic structure in the leads and in the device re- 
gion typically based on the ground state self-consistent 
Kohn-Sham orbitals. Because of the conceptual appeal of 
this approach, the earlier works were not focused on its 
formal justification. Recently, however, several studies 
have provided the rigorous theoretical foundation for the 
above DFT-based methodologies calculating the steady 
currents in open electronic system s 34 ' 43 ! 44 ' 45 . Therefore, 
the question concerning the validity of the present ap- 
proach relies mostly on the proper description of the ex- 
change and correlation within the DFT approximation. 
We focus below only on two aspects of the DFT approach 
that seem to be the most important for understanding of 
the discrepancy between the calculations and the exper- 
iment. 

The exchange and correlations are commonly 
accounted for within the local spin density 
appromaximatioiiS 3 .. As mentioned in preceding 
sections, for two dimensional systems there arc two 
commonly used parameterizations, namely the param- 
eterization of Tanatar and Ceperley 2 ^ and Attacalite 
et al2^. The validity of these approximations was 
tested for few-electron quantum dot systems and, 
generally, a very good agreement agreement with the 
exact diagonalization and/or variational Monte Carlo 
calculations is found 4 ^. Taking into account that these 
parametrizations give practically the same results for the 
QPC conductance, we do not expect the utilization of 
the above parametrizations to be a source of a significant 
discrepancy between the calculated conductance and the 
experiment. Instead, we focus on another aspect of the 
choice of the exchange energy functional which arises in 
the systems with a variable particle number. (Note that 
the QPC structure, being an essentially open structure, 
belongs to this class of systems). This aspect is related 
to the infamous "derivative discontinuity problem" of 
the DFT originating from the discontinuous dependence 
of V xc on the particle number 23 . (Note that the LDA 
does not include any derivative discontinuity in the V xc ). 

The problem with the derivative discontinuity (lead- 
ing to uncompensated self-interaction corrections in the 
V xc ) has been recently recognized in the standard quan- 
tum mechanical transport calculations in molecular sys- 
tems and atomic wire o 43 ' 47 ' 48 . Typically, such the calcu- 
lations provide a good quantitative agreement with the 
experimental data for systems where the coupling to the 
leads is strong and the conductance exceeds the conduc- 
tance unit Go = 2e 2 /h (for example, for atomistic metal- 
lic wires and related systems). At the same time, for 
weakly coupled systems such as organic molecules the 
standard DFT+LDA approach leads to the orders-of- 
magnitude discrepancy between the measured and cal- 
culated currents and to incorrect predictions of the con- 
ducting (instead of experimentally observed insulated) 



phas o ' ' ' ' 50 . There have been attempts to explain 
these discrepancies by insufficient modelling, such as 
atomistic structures or contact coupling. Several recent 
studies however, attributed this discrepancy to more fun- 
damental reasons, identifying the lack of the derivative 
discontinuity in LDA as a major source of error in the 
DFT-based transport calculation s 47 ' 48 . For example, To- 
her et al£l argue that LDA approximation is not suit- 
able for transport calculations for the case of weak cou- 
pling. They propose a simple corrective scheme based on 
the removal of the atomic self-interaction. Their correc- 
tive scheme restores an agrement with the experiment, 
opening a conduction gap of the I- V characteristics of a 
molecular junction instead of a metallic behavior follow- 
ing from the standard DFT+LDA approach. (It is in- 
teresting to note that similar self-interaction corrections 
practically do not effect the conductance and electronic 
density for the case of strongly-coupled systems 4 - 7 ). Al- 
ternative corrective schemes replacing the Kohn-Sham 
data for the weakly coupled region with their counter- 
parts obtained from a Hartree-Fock analysis (taking care 
of the self-interaction problem) are suggested by Evers et 
al££ and Palacios^ 9 -. Note that various approaches to the 
description of quantum transport for the case of weakly 
coupled systems (accounting for the charge quantization 
and thus eliminating the self-interaction errors) were dis- 
cussed in Refs. 50.51.j52l 

We argue here that a similar problem related to the 
derivative discontinuity may be the reason why the stan- 
dard DFT+LDA approach fail to describe the observed 
0.7-anomaly in the QPC. Indeed, the formation of a spin- 
polarized charge droplet predicted within the DFT+LDA 
approach implies that electrons are trapped in weakly 
coupled quasi-bound states in the center of the QPC. As 
mentioned above, in the case of weak coupling the lack 
of the derivative discontinuity in the LDA approximation 
causes the orders-of-magnitude discrepancies between the 
theory and experiment for the molecular systems. Thus, 
it would be reasonable to expect that due to the same 
reason the LDA approximation is not suitable for the 
case of the QPC structure as well. Because the correc- 
tive schemes accounting for the derivative discontinuity 
are shown to strongly affect the electron density and the 
energy levels in the system, and because of the apparent 
failure of the standard DFT+LDA approach to repro- 
duce the 0.7 anomaly, we conclude that the formation of 
the magnetic impurities in the QPC might be an artefact 
of the LDA due to the lack of the derivative discontinuity 
related to the spurious self-interaction. 

Based on the above discussion we conclude that an ac- 
curate description of the QPC structure might require 
approaches that go beyond the standard DFT+LDA 
scheme and account for the derivative discontinuity in 
the V xc or utilize similar corrective schemes eliminat- 
ing the self-interaction errors of the DFT+LDA. It is 
not clear at the moment whether the recipes for the ac- 
counting of the derivative discontinuity and eliminating 
the self-interaction correction developed for the molecu- 
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lar junction o i ' can be adapted for the QPC struc- 
ture. Such the corrective schemes are remained to be 
implemented and it remains to be seen whether this can 
bring the calculated conductance to the closer agreement 
with the experiment. 

As an indirect support of the above arguments we 
notice that similar spin DFT conductance calculations 
(within the same LDA approach and the same param- 
eterization of V xc ) 53 reproduce quantitatively the mea- 
sured spin-resolved magnetoconductance of the QPCs in 
the integer quantum Hall regime 54 . In this case the edge 
state regime is reached such that the transport through 
the QPC correspond to the strong coupling regime. 

Finally, a comment is in order concerning the Kondo 
physics. The Kondo effect was suggested by Cronen- 
wett et a/£ as a possible source of the 0.7 anomaly; 
at the same time, the experimental studies of Graham 9 
Rokhinson 1 ^ seem to rule out this interpretation. We 
stress that the present mean-field approach based on the 
standard DFT+LDA formulation is not able to address 
this effect. Moreover, we do not believe that the pre- 
diction of the magnetic impurity formation within the 
DFT+LDA approach can support or rule out the Kondo 
physics in the QPC. Indeed, the Kondo enhanced con- 
ductance channels (if they exist) would change the elec- 
tron density inside the QPC (beyond that one predicted 
by the DFT). Thus, the DFT+LDA predictions for the 
equilibrium electron density in the QPC should be cor- 
rected in a self-consistent way to account for this addi- 
tional density, and it is not obvious whether the spin- 
polarized quasi-bound states predicted in the framework 
of the DFT+LDA would survive this correction. 



VI. CONCLUSION 

We have developed an approach for full quantum me- 
chanical transport conductance calculations in open lat- 
eral split-gate structures that starts from the lithographi- 
cal layout of the device and is free from phenomenological 
parameters like coupling strengths, charging constants 
etc. The whole structure, including the semi-infinitive 
leads, is treated on the same footing (i.e. the electron- 
electron interaction is accounted for in both the leads 
as well as in the central device region). The electron- 
electron interaction and spin effects are included within 
the spin density functional theory in the local spin den- 
sity approximation, and the conductance calculated on 
the basis of the standard Green's function technique. 

The developed method was applied to calculate the 
spin-resolved conductance through a QPC. Close to the 
pinch off a spin-degeneracy of the spin-up and spin-down 



conductance channels is lifted and the total conductance 
shows a broad feature peaked at 0.5 x 2e 2 /h (correspond- 
ing to a complete suppression of one of the spin chan- 
nels). A similar feature is also present in the range of 
the gate voltages where a second step in the conduc- 
tance develops. The lifting of the spin-degeneracy and 
the suppression of one of the spin channels are the generic 
features of all studied QPCs (both very short and very 
long; 40 < I < 500 ran). The calculated conductance also 
shows a hysteresis for forward- and backward sweeps of 
the magnetic fields. These features in the conductance 
are related to the formation of weakly coupled (quasi- 
bound) states in the constriction of the QPC (predicted 
in previous DFT-based studie a 16 i 17 ) . A comparison of the 
obtained results with the experimental data shows how- 
ever, that while the spin-DFT based "first-principle" cal- 
culations predict the spin polarization in the QPC struc- 
ture, the calculated conductance clearly does not repro- 
duce the 0. 7 anomaly observed in almost all QPCs of 
various geometries. 

In order to understand why the calculated conductance 
fails to reproduce the 0.7 anomaly, we critically exam- 
ine the major features of the DFT-based conductance 
calculations. We suggest that inability of the standard 
DFT+LDA approximation to reproduce the 0.7 anomaly 
might be related to infamous derivative discontinuity 
problem of the DFT leading to spurious self-interaction 
errors not corrected in the standard LDA 23 . This prob- 
lem has been recently recognized in similar DFT-based 
transport calculations for the molecular junctions (show- 
ing orders-of-magnitude discrepancies with the experi- 
ment) where it has been demonstrated that the LDA ap- 
proximation is not suitable for transport calculations for 
the case of weak couplin g 47 ! 48 . We thus conclude that 
the formation of the weakly coupled magnetic impurities 
in the QPC might be an artefact of the LDA due to the 
lack of the derivative discontinuity and related spurious 
self-interaction. Our results suggest that an accurate de- 
scription of the QPC structure requires approaches that 
go beyond the standard DFT+LDA scheme and that ac- 
count for the derivative discontinuity in the V xc or utilize 
similar corrective schemes eliminating the self-interaction 
errors of the DFT+LDA. 
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